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Abstract: We present a new parton-shower algorithm. Borrowing from the basic ideas of 
dipole cascades, the evolution variable is judiciously chosen as the transverse 
momentum in the soft limit. This leads to a very simple analytic structure of the 
evolution. A weighting algorithm is implemented, that allows to consistently treat 
potentially negative values of the splitting functions and the parton distributions. 
We provide two independent, publicly available implementations for the two event 
generators PYTHIA and SHERPA. 


1 Introduction 

Parton showers and fragmentation models have been used for more than three decades to predict the dy¬ 
namics of multi-particle final states in collider experiments 0. More recently, the traditional approaches 
implemented in HERWIG [2], PYTHIA |TT| and SHERPA [3] were supplemented by methods based on dipole and 
antenna factorization [311130 E]- A characteristic feature of these new shower programs is the description 
of QCD coherence in the color dipole picture [10], which has first been implemented in the ARIADNE Monte 
Carlo m- In this article we present a dipole-like parton shower similar to existing ones, but we focus on 
the simplest implementation and enforce sum rules and DGLAP collinear anomalous dimensions. We choose 
ordering variables based on transverse momenta in the soft approximation, while existing dipole-like shower 
models employ collinear transverse momenta. As such, the model is a hybrid of dipole and parton shower. 
These choices will eventually allow to compare with analytic approaches, such as CSS [12] and SCET [l3l . 
In the past decade, the matching of parton showers to NLO calculations mug and the merging of LO [IB] 
and NLO matched results for different jet multiplicity Hang was in the focus of interest of the majority of 
Monte-Carlo developers m- Comparably few efforts were made to provide publicly available implementa¬ 
tions of parton-showers 0000 or to improve their formal accuracy [20] , and even fewer of the new parton 
showers have made their way into complete event generators used by experiments. When comparing results 
of matched and merged calculations, it is therefore often unclear whether a particular difference stems from 
mismodeling in the parton shower, from differences in the matching or merging algorithm, or simply from 
technical problems. Similarly, when comparing the results of different event generators at the hadron level 
it is often unclear whether differences should be ascribed to the hadronization model, to the simulation of 
multiple scattering / rescattering effects, or to the parton shower. We intend to remedy this situation to some 
extent, by providing two implementations of one and the same algorithm, to be used with the two different 
event generation frameworks PYTHIA m and SHERPA [22]. We subject our codes to rigorous scrutiny by 
comparing their predictions at the sub-permille level. 

This paper is organized as follows: Section [2] reviews the basic parton-shower formalism. Section [2] explains 
the construction principles of our new parton shower, which we call DIRE (acronym for Dipole REsummation). 
Section [4] contains the validation of the numerical implementation, and Sec. 0 presents a comparison of the 
predictions from DIRE with experimental measurements. Section 0 contains some concluding remarks. 


2 Parton-shower formalism 


The evolution of parton densities and fragmentation functions in the collinear limit is governed by the 
DGLAP equations [25] : 
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where P ab are the regularized evolution kernels. Assume that we define P a b in terms of unregularized kernels, 
P a b, restricted to all but an e-environment around the soft-collinear pole, plus an endpoint contribution. 
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For finite e, the endpoint subtraction can be interpreted as the approximate virtual plus unresolved real 
corrections, which are included in the parton shower by enforcing unitarity. The precise value of e is defined 
in terms of an infrared cutoff on the evolution variable, using four-momentum conservation. When ignoring 
momentum conservation, this cutoff can be taken to zero, which allows us to identify [P ba {z)] + as the e —> 0 
limit of P ba (z,e). For 0 < e <C 1, Eq. (2.11 changes to 
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one can define the generating functional for splittings of parton a as 
Pa(x,t,l-l 2 ) = fa{x, t)A a (t, /i 2 ) = fa(x,ig 2 )n a (x,t,H 2 ) , 
where 
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In this context, n a (:r, t, /i 2 ) is the probability that the parton does not undergo a branching process between 
the two scales /r 2 and t [23]. Equation (2.3) can now be written in the simple form 
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The generalization to an n-parton state can involve multiple PDFs and fragmentation functions: 
dhi P s ($„,t,ii 2 ) 
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In this context, we have extended the argument of the generating functional to <F n , which denotes the n- 
parton phase-space configuration, including all light-cone momentum fractions, x-i and Zj. for initial-state 
(IS) and final-state (FS) partons. T also depends on all parton flavors, denoted by a. If we do not fix 
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Figure 1: Kinematics in the initial-state parton splitting process a —> {aj}j. 


the momenta of final-state hadrons, the fragmentation functions can be integrated over Zj, leading to the 
simplified formula 
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The change from $ n to 4> n signals that F has become independent of Zj. An observable-dependent generating 
functional for the parton shower can now be defined recursively as 


f 4 ' dt 
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dint 
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( 2 . 10 ) 


The first term includes all virtual corrections and unresolved real emissions, resummed into a no-branching 
probability. The second term describes a single branching, followed by further parton evolution. Both 
terms can be generated simultaneously by implementing the veto algorithm j25j for Eq. (2.6). We have 


introduced an observable, O, that measures the kinematics of the final state. In general, this observable will 


act differently on the no-emission term and on the emission term. In the trivial case that 0 = 1, Eq. (2.10) 
returns the unitarity constraint, > n ,t,t'; 1) = 1. 

Generating a branching in the parton shower involves selecting a new color topology for the n + 1-particle 
state. For non-trivial color configurations, T will therefore depend on the color assignment in the large-7V c 
limit. While it is in principle necessary to keep track of this dependence, we omit any notation relating to 
color in order to simplify our final formulae. The selection of color topologies proceeds as in existing dipole¬ 
like parton showers, which is described in great detail in |Bj. It is straightforward to extend our notation in 
this regard. 

The choice of evolution variable is crucial. At leading color the soft radiation pattern emerges from the 
coherent gluon radiation off “color dipoles” that are spanned by the two partons at opposite ends of a color 
string |2B] ■ This mandates the choice of an evolution variable which treats these two partons democratically. 
In other words, the evolution variable should be identical no matter whether one or the other of two color- 
connected partons is considered the radiator. This will be the guiding principle for its selection in Sec. [3] 
The splitting functions, P(z), are formally defined in the collinear limit, and they do not reflect the soft 
radiation pattern outside the collinear region. Traditionally, this problem is dealt with by imposing angular 
ordering constraints on the final-state phase space (5]. Alternatively one can use the approach of Catani 
and Seymour [27], and introduce a f-dependence in the splitting functions that restores the correct soft 
anomalous dimension at one-loop order ii- We will use this procedure in the next section. It is important 
that the modified splitting functions satisfy the sum rules, which are enforced by Eq. (2.2) and by the 
corresponding flavor sum rule [24] , The new splitting functions may also be negative in the non-singular 
phase-space region. This requires a modification of the Sudakov veto algorithm [MUMUm, and it entails 
an analytic event weight to ensure that both emission- and no-emission probabilities are accounted for. We 
find that, in our parton-shower approach, the variance of this weight is small. In fact, for both final-final 
and initial-initial dipoles momentum conservation guarantees that no negative weights can arise from the 
splitting functions. Negative weights may however appear in initial-initial configurations due to negative 
values of the PDFs. 
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3 Construction of the Dire shower 


A basic branching process is sketched in Fig. [T] In this case we consider initial-state evolution. We employ 
the kinematics from Ref. 2~7 l.’llll, which we summarize in Appendix |A| For initial-state splitters with initial- 
state spectator, all particles typically have zero on-shell mass, which greatly simplifies the calculation. Their 
momenta can be parametrized in terms of the light-cone momenta p a and pb , using the standard Sudakov 
decomposition m 
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In the Catani-Seymour approach m, the correct soft anomalous dimension is obtained after replacing the 
soft enhanced term in the splitting functions by a partial fraction of the soft eikonal for the color dipole 
defined by the splitting parton and its spectator. Schematically this can be done as follows: 
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If we define the evolution variable of our parton shower to be a scaled transverse momentum, t = (z — v) k\, 
then the soft-enhanced term in Eq. (3.2) is conveniently expressed as 
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Note that the evolution variable has the desired symmetry property, i.e. it is symmetric in emitter and 
spectator momentum. More precisely, our evolution variable is the exact inverse of the soft eikonal. As 
such, it is different from the hardness parameter, k\. Consequently, the parton shower will fill the entire 
final-state phase space, even for factorization scales much smaller than the hadronic center of mass energy. 


We define the evolution kernels using the replacement of the soft enhanced term in Eq. (3.2). Additionally, 


we require the collinear anomalous dimension to be unchanged. Imposing flavor and momentum sum rules, 
these two requirements determine the complete set of leading-order spin-averaged splitting functions as: 
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P qq (z,n 2 ) = 2C F 

P gg (z,n 2 ) = 2C A 

P qg (z,n 2 ) = 2C F 

P gq (z, K 2 )=T R z 2 + (l-z) 2 

The corresponding anomalous dimensions are listed in Appendix [Bj Using the phase-space factorization 
derived in EH, we obtain the following differential branching probability: 
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where 2 z± = 1 + x ± \J{1 — x) 2 — 4k 2 Note that this implies x/(z — v) < 1, i.e. the light-cone momentum 
fraction entering the PDFs is well defined. This variable has changed compared to Eq. (2.1) to account 
for four-momentum conservation, while Eq. (2.1) remains valid in the collinear limit, v —> 0. In addition, 


the strong coupling is evaluated at the evolution variable, hence the Landau pole is avoided by the infrared 
cutoff of the parton shower, to (which is of order 1 GeV). Finally, the splitting kinematics are constructed 


as described in App. A.5 


The technical implementation of Eqs. (3.4) in terms of dipole terms proceeds as described in [B]. That is, 
we divide the splitting function according to the number of spectator partons in the large-A(. limit, and 


1 The actual value of the integration limits on z does not have to be computed explicitly. In practice, one generates 

Monte-Carlo events in the maximum range x < z < 1, and vetoes events violating momentum conservation (cf. [25]). 


4 



































sum over color-adjacent splitter-spectator pairs. The corresponding evolution equations are straightforward 
extensions of Eq. (3.51, and therefore we do not list them here. The same reasoning applies to all dipole 
types discussed in the following. 

Initial-state splitters with final-state spectator are treated along the same lines. The construction of final- 
state momenta is described in App. |A.3[ The kinematics are now defined in terms of the following variables 
(pk is the final-state spectator momentum) 


2 = 1 - 


PjPk 


u = 


PjPa 


PaPj + PaPk PaPj + PaPk 

We define the evolution variable, t, and its dimensionless variant, k 2 , as 
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The evolution variable is symmetric in emitter and spectator momentum, up to a rescaling by I/ 2 2 . The 
differential branching probability is: 
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Final-state splitters with initial-state spectator must have the same evolution variable due to symmetry 
arguments. Therefore the asymmetric rescaling by I/ 2 2 in the IF case must also be applied in the FI case. 
The only choice to be made concerns the splitting variable, which is taken to be 1 — u. The differential 
branching probability is: 
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where 2 = 1 — t/Q 2 /( 1 — 2 ). The integration limits are given by 2 _ = 1 — t/Q 2 /( 1 — x) and z + = 1 — t/Q 2 . 
The construction of final-state momenta is described in App. |A.2| 

For final-state splittings, ij —>• i,j, with final-state spectator, k, we use the variables 
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The symmetric evolution variable, its scaled variant k 2 and the splitting variable are defined as 
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The differential branching probability is: 
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where 2 z± = 1 ± \J\ — 4t 0 /Q 2 . The splitting kinematics are described in App. 

Note that the scaled transverse momentum defined in Eq. (3.11) is substantially different from the ones 
defined in u®, which can be written as 


~2 _ _ ZpiPj PiPkPkPj 

Q 2 Q 2 ({Pi + Pj)pk ) 2 ' 


(3.13) 


This variable is symmetric in i and j, but not in i and k , which would be required in order to interpret it 
as the inverse soft eikonal for gluon radiation off the dipole spanned by i and k. Kinematically, Eq. (3.13) 
represents the transverse momentum of partons i / j with respect to the anti-collinear direction defined 
by k. This is what we call a “collinear” transverse momentum. In contrast, Eq. (3.11) can be interpreted 
as the transverse momentum of the two daughter dipoles (ij) / ( kj) in the center-of-mass frame of the 
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decaying dipole m- In this case, i defines the collinear, and k defines the anti-collinear direction, making 
the symmetry explicit. We refer to such a definition as a “soft” transverse momentum. 

The change in the definition of transverse momenta compared to existing py-ordered dipole-like parton 
showers mi also involves changing the splitting variable, in order to reduce the related Jacobians to unity 
while maintaining Eqs. (3.41, simultaneously for all dipole types. In contrast, the kinematics mapping is 


identical to the previously published methods @ 22 ]. 

If massive quarks are involved in the branching process, we would like to map the evolution variable to the 
soft enhanced term of the full matrix element, just like in the massless limit. The singularity structure in 
the soft limit is given in m■ For the most involved case of two massive radiators, i and k, it leads to an 
eikonal of the form 


PiPk 


2 (PkPj) 2 
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The inverse of this function is difficult to interpret. Its scaling property in the soft limit, however, is 


completely determined by the first term in Eq. (3.14), whose inverse can therefore be used to define an 


ordering variable for the evolution of massive partons 
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with 7 defined in Appendix |Aj and with s lk = (pi + p k ) 2 - Here we have defined the massless equivalent of 
the evolution variable, k^ 2 , and the light-cone momentum fractions, ^ and Q in a Sudakov decomposition 
of the gluon momentum, pj , along the directions of pi and p k : 
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Equation ( 3.15| ) is valid in the soft limit. For practical purposes the denominator ptpk in the evolution 
variable should be the hard scale of the radiating dipole, which is given by (2 p a Pb) 2 /Q 2 i (ZPaiPi +Pk)) 2 /Q 2 
and Q 2 for II, IF/FI and FF dipoles, respectively. 


The splitting functions for massive partons can be taken from Eq. (3.4) and be modified according to 


We use the following unregularized massive kernels for final-state splitter with final- or initial-state spectator: 
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The relative velocity between two momenta, p and q , is defined as 
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and v~ j. and v-j 3 _k stand for the relative velocities between the emitter parton and the spectator before and af¬ 
ter the branching, respectively. The branching probabilities are modified as P a b(z, k 2 ) —> J(z , n 2 ) P a b(z , k 2 ), 
where J(z,t) is a spectator-dependent Jacobian factor |30([6]. It is unity, except for the case of final-state 
splitter with final-state spectator, where 
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using Q 2 = 2 pp k + 2 (pi +Pk)Pj, cf. Eq. (3.10). The phase-space boundaries are given by the roots of the 
Gram determinant 


4A 3 = 2 ptpj 2 pjp k 2 pip k - (2 pipj) 2 ml - (2 pjp k ) 2 m 2 - (2 PiPk) 2 m 2 + 4 m 2 m 2 m 2 k . 


(3.20) 
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While the massless case leads to simple constraints on z, the general massive case generates a rather involved 
functional form of the ^-boundary as a function of t. Algorithmically, it is preferable to use the veto 


algorithm [25] to implement this constraint, or to use Eqs. (3.11) and (3.7) and evaluate the constraint in 
collinear variables, where m 


Z± = 


Pipj 


■ m,; 


(pi+PjY 


(i± 


V ij,j V ij,k , 


(3.21) 


In final-state splittings with initial-state spectator the PDF is evaluated at x/z/{ 1 + (m 2 j — m 2 — m 2 )/Q 2 ). 


Correspondingly, the theta function in Eq. (3.9) changes to Q(Q 2 (1 — z)(l — x Q 2 /(Q 2 + m 2 j — m 2 — m 2 )) — t). 
For initial-state splitter with final-state spectator the mass-dependent splitting functions are 


PU(Z,K 2 ) = C F 
pV{z,n 2 ) = 2C A 


Z z + K z Q z 1 — u 



(3.22) 


4 Validation 

In this section we validate the numerical implementation of the DIRE parton shower. The two event generation 
frameworks PYTHIA ^3] and SHERPA 0] are used to construct two entirely independent Monte Carlo programs. 
Aside from a thorough cross-check of the implementation, this allows, for the first time, to simulate Deep 
Inelastic Scattering in PYTHIA 8. We employ the CTlOnlo PDF set [33], and the corresponding value of the 
strong coupling. Following standard practice to improve the logarithmic accuracy of the parton shower, the 
soft enhanced term of the splitting functions is rescaled by l + a s (t)/(2n)K, where K = (67/18 —7 t 2 /6) C a — 
10/9 T R n f [31]. 

Figures and [4] each show a comparison between the results from DlRE+PYTHIA and DlRE+SHERPA. 
Each simulation contains 10 8 events. The lower panels present the deviation between the two predictions, 
normalized to the statistical uncertainty of DlRE+SHERPA in the respective bin. This distribution should 
exhibit statistical fluctuations only. We validate quark splitting functions in the reactions e + e _ —> hadrons 
(Fig. | left), e + p —¥ e + jet (Fig. [3] left), and pp —> e + e (Fig. [4] left). These three cases cover all possible 
dipole types, i.e. final-state splitter with final-state spectator, final-state splitter with initial-state spectator 
and vice versa, and initial-state splitter with initial-state spectator. Gluon splitting functions are validated 
in the reactions t + t - —>• hadrons (Fig. [ 2 ] right), r + p —► r + jet (Fig. [ 3 ] right), and pp —> t + t~ (Fig. [ 4 ] left), 
all mediated by Higgs-boson exchange using HEFT [35j . 
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Figure 2: Validation in e + e —► hadrons and t + t 
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Differential kj -jet resolution at parton level (Breit frame) 
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Figure 3: Validation in e q —> e <7 and t g ^ t g 
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Figure 4: Validation in gg — > e + e and gg — > t + t 
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5 Results 


In this section we compare DIRE predictions from SHERPA [22] with experimental data. When applicable, 
we use the CTlOnlo PDF set [33] and the corresponding strong coupling. We employ the kinematics scheme 
from App. |A.3| Our results include the simulation of QED radiation in the case of Drcll-Yan lepton pair 
production [35], and hadronization in the case of e + e _ —>• hadrons 133 ■ Otherwise they are given at the 
parton level in order to exhibit the features of the DIRE shower only. Analyses are performed with RIVET J3B] . 

Figure [5] shows predictions from the DIRE parton shower for differential jet rates in the Durham scheme 
compared to experimental results from the JADE and OPAL collaborations £351- The perturbative region is 
to the right of the plot, and y ~ 2.8-10 -3 corresponds to the &-quark mass. The simulation of nonperturbative 
effects dominates the predictions below ~ 10~ 4 . We observe that, in the perturbative region, the results are 
in excellent agreement with the experimental measurements. 

Figure [6] shows a comparison for event shapes measured by the ALEPH collaboration [ID] . The perturbative 
region is to the right of the plot, except for the thrust distribution, where it is to the left. We notice some 
deviation in the predictions for jet broadening and for the C-parameter. However, these deviations are mostly 
within the 2a uncertainty of the experimental measurements, and they occur close to the nonperturbative 
region. It can also be expected that the simulations improve upon including matrix-element corrections 
or when merging the DIRE shower with higher-multiplicity calculations. This has been demonstrated, for 
example, in El- 

Figure [7] shows angular correlations in comparison to ATLAS data from EH, and the transverse momentum 
spectrum of the Drell-Yan lepton pair in comparison to ATLAS data from [42] ■ It is well known that 
pure parton shower predictions are insufficient to describe these measurements. Therefore, we merge our 
parton shower with 1-jet matrix elements using the CKKW-L procedure M- In order to assess the related 
uncertainties, we vary the merging cut by a factor 2 around the central value of Q C ut = 10 GeV. The 
associated uncertainty band is shown in light red. The size of the variation is comparable to the statistical 
uncertainties, which are displayed as error bars on the Monte-Carlo prediction. 

Figure [8] shows di-jet azimuthal decorrelations in different regions of jet transverse momentum. We compare 
DIRE predictions with experimental results from the ATLAS collaboration j33]. This observable tests for 
higher-order effects in some detail [44] . 


6 Conclusions 

We presented a new dipole-like parton-shower algorithm, that is constructed along very simple arguments: 
Firstly, the ordering variable should exhibit a symmetry in emitter and spectator momenta, such that the 
dipole-like picture can be re-interpreted as a dipole picture in the soft limit. At the same time, the splitting 
functions are regularized in the soft anti-collinear region using partial fractioning of the soft eikonal in the 
Catani-Seymour approach. They are then modified to satisfy the ordinary sum rules in the collinear limit. 
This leads to an invariant formulation of the parton-shower algorithm, which is in complete analogy to the 
standard DGLAP case. We computed the anomalous dimensions, which match previous results for dipole¬ 
like parton showers. We presented first phenomenologically relevant predictions using the new algorithm, 
and we observe very good agreement with experimental data from LEP and LHC experiments. 
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Figure 5: DIRE predictions in comparison to LEP data from [39]. 
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Figure 8: DIRE predictions in comparison to ATLAS data from |4.‘j| . 


14 



































































































A Parton-shower kinematics 


The precise algorithm for constructing the splitting kinematics depends on the type of splitter and spectator 
parton. There are four separate cases. Note that initial-state partons are assumed to be massless for 
collinear PDF evolution to be valid. In practically implemented parton-shower algorithms they are often 
taken massive instead, in order to obtain a better description of experimental data. Therefore we give the 
kinematics formulae with full mass-dependence, including initial-state parton massesj^] 


A.l Final-state splitter with final-state spectator 


1. Determine the new momentum of the spectator parton as 


u 

Pk = 


Pk 


‘mv 


A (q 2 ,s ij ,m 2 k ) 


A (q 2 , , to 


•j? 

n' 


2 ) 


g 2 + m\ - Sjj 
2 q 2 q 


with A denoting the Kallen function A (a, 6, c) = (a — b — c) 2 —4 be 
and Sij = y ij}k ( q 2 - m 2 k ) + (1 - y ijtk ) (m 2 + m 2 ). 

2. Construct the new momentum of the emitter parton, p,, as 

m = _ 7(g 2 ,%,Wfc)<- - Sij pi to 2 + k 2 ± p^-m 2 /'r{g 2 ,s lJ ,m 2 )p^ 

/3(q 2 ,Sij,ml) Zi /3{q 2 ,s ij ,ml) '- L 


where q = p tj + p k , (3(a, b , c) = sgn (a - b - c) y/\(a,b,c), 

2 7 ( 0 , b,c) = (a- b - c) + /3(a , 6 , c) and pC = q^ - p k . 

The parameters z* and kj_ = — k 2 ^_ of this decomposition are given by 


T - s ij - m l 


P (g 2 P 


+m 2 — m 2 


7(g 2 ,- 


,ml) q 2 - - 


kl 


= Zi (1 - Zi) Sij - (1 - Zi) m 2 i - 


Zi ru¬ 


in the massless case, this algorithm reduces to the following [ 6 ] 

Pi = ZiPij + Vij,k{ 1 - Zi)Pk + K , 

Pj = (! - Zi)Pij + Vij,kZiP k - k£ , 

Pk = { l ~yij,k)Pk » 

where Z; = (z - Vij,k)/{ 1 - Vij,k) and y ijtk = k 2 /( 1 - z), cf. Sec. [s] 


(A.l) 


(A.2) 


(A.3) 


(A.4) 


A.2 Final-state splitter with initial-state spectator 

1. Determine the new momentum of the spectator parton as 


P 


m _ 

a 



q-Pa A I A(g 2 ,s lj ,m 2 a )^4m 2 a p 2 J± 
q f\ g " / V A (g 2 >TO 2 i ,m2) -4m 2 a p 2 J± 


q 2 + m 2 a - Sij 


lq \\ 


q \\ 


(A.5) 


where q = p a - Pij, q\\=q + p ij± and = (1 - 1 /x ijt a) (q 2 - m 2 a ) + (m? + m?)/xi j}a . 

2. Proceed as in Sec. B except that m k —> m a and p k —>■ —p a - 

2 We use the notation and definitions of im both in the massive and in the massless case. This implies z, = {pip k )/(PiPk + 
PjPk ) and yij ik = 2piPj/Q 2 for final-final dipoles, Zi = (piPa)/{PiPa + PjPa) and 2 Xij, a = Q 2 /(piPa + PjPa) for final-initial 
dipoles, Uj = PjPa/{PjPa + PkPa) and 2xj k , a = Q 2 /{PjPa + PkPa) for initial-final dipoles, and Xj^ a b = Q 2 / (?p a Pb) and 
Vj = (p a Pj)/{.PaPb) for initial-initial dipoles. 
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In the massless case, this algorithm reduces to the following [B] 
1 


Pi = ZiPij + (i - h) 




_L ’ 


°ij,a 


Pi =(1- Zi) Pii + Zi - -^ ~ kX 


_L ’ 


v ij,a 


p£ = 


■Pa , 


v ij,a 


where Zi = z and Xi ha = 1 — k 2 /(1 — 2 ), cf. Sec. [ 3 ] 


(A.6) 


A.3 Initial-state splitter with final-state spectator (local recoil) 

1. Determine the new momentum of the splitting parton as 


Pa = Paj - 


<1-Pa3 

.2 9 || 


\{q 2 ,s jk ,m 2 a ) -Am 2 a p 2 kl _ q 2 +m 2 a -Sj k 


n 


A(g 2 ,m>.h)_4m2 p2 ± 


' q \\ 


% 


where q = p aj - p k , q^=q + p k± and s jk = (1 - l/a^ fc)0 ) (g 2 - m 2 ) + (?n 2 + m\)/x ikta . 


2. Proceed as in Sec. E3 except that Zj —)• Uj , m k —► m. 

In the massless case, this algorithm reduces to the following [6, 


(A.7) 


rrij 


m k , p k —» ~Pa and pj -> p k . 


p£ = 


%jk,a 


a 3 


1 %j k , 


a 


Pj = (! - u j) „ ' ' Paj 
1 

/i _ f — %jk,a ~fi 




_L 5 


% = % 


%jk,a 

where Xjk, a = 2 and w 


%jk, 

Paj + i 1 ~ U j) Pk + k ± > 


(A.8) 


= k 2 /(1 — 2 :), cf. Sec. [ 3 ] 


A.4 Initial-state splitter with final-state spectator (global recoil) 

This scheme can be chosen as an alternative to the one in Sec. |A.3| The algorithm is equivalent to the method 
outlined in [281I32j . It alleviates a formal problem with transverse momentum resummation in Drell-Yan 
type processes [IS], but it is numerically less stable, and therefore employed only if | Xj ka — uj\ > e, where 
e ~ 1CT 4 . 


1. Determine the new momentum of the spectator parton as 


Pk = Pk 


1 ’Pk 


I A(q 2 ,s aj -,m%) 
A (q 2 ,mlj,ml) 


g 2 + m 2 k -s a j 
2 q 2 q 


where q = p k - p aj and s aj = Uj/xj k , a ( q 2 - m 2 k ) + (m 2 + m 2 ) (1 - Uj)/x jk , a . 
2. Construct the momentum of the emitted particle, pj , as 


Pj = - Z j 


_ l{q 2 ,Saj,rnl)Pg 1 + Sajpg rrij + k 2 p k + m k / r y(q 2 , s a j,m k ) p t 


P{q 2 


,Saj,m 2 k 


/3(q 2 ,' 


•aj, "^ k ) 


The parameters Zj and = — k\ of this decomposition are given by 


Zi = 


q 2 - Saj - m k Xj k 5 a - 1 


P(q 2 ,s 2 ai ,m 2 ) 


ml 


Saj + TO 2 — m 2 


Xjk,a - Ui 7 (q 2 , s aj , ml) q 2 - Saj - m 2 k 


k i = Z j (! - Zj) Saj - (1 - Zj) TO 2 - Zj m 2 a 


K, 


(A.9) 


(A.10) 


(A.ll) 
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3. Boost p a and all final state particles into the frame where p a is aligned along the beam direction, 
with p b , the opposite-side beam particle, unchanged. 

In the massless case, this algorithm reduces to the following 


a _ 1 ~ u 3 


Pa = 


Pi = 


Xjk,a 
1 %jk,a 

Xjk,a 


V ■ + 
Vclj ' 


1 Xjk,a 


Xjk ,a Xjk ,a u*j 


Pk + k X . 


zr. + 
Vclj ' 


'Jjj 1 Uj u 

- -— Pk + k ±> 

Xij,a Xjk, a - Uj 


(A.12) 


-;£)*• 


Note that, following the arguments in [9], the light-cone momentum fraction entering the PDF is still given 

by x/z = x/xij'k- 

A.5 Initial-state splitter with initial-state spectator 

1. Determine the new momentum of the splitting parton as 


pH = pH- 


7 



^(^abi 






7 (Sab,m 2 a ,m 2 b ) 


pH 


(A.13) 


where q = p a + Pb and s ab = ( q 2 - mf)/xj tab + (1 - l/x^ ab ) (ml + ml). 

2. Construct the momentum of the emitted parton, pj, as 

7 (s a b,m 2 a ,ml)p£ -m 2 aP H m 2 + k 2 ± p g - ml/^(s ab , m 2 , m 2 b )p£ 


pH = (1 ~ z a 


P(Sab, i 


m b) 


1 ~ Z, 


aj j3(s ab ,m 2 a ,ml 


-kH 


_L ! 


The parameters z a j and = — k\ of this decomposition are given by 


(A.14) 


Sab - ml- ml 


°aj 


/3(s ab ,m*,m 

Zaj 


l,ml) 

aj 


(Xj,ab + v j ) 


S aj + ml — m^ 


k — Z a j (1 Zaj) m a ( X Zaj) Saj Z a i m 


7 (Sab, ml, ml) s ab -m 2 a -m 
aj ) J aj — Z a j VTlj , 


(A.15) 


3. Boost all remaining final-state particles into the frame defined by p a + p b — Pj, using the algorithm 
defined in Sec. (5.5) of [27]. 

In the massless case, this algorithm reduces to the following [6. 

1 


Pa = 


~LL 

Paj 


Xj^ab 

1 Xa nh Va 


Pj = 


Paj+VjP b ~K , 


(A.16) 


Xj,ab 

fcf = A (p ai + p b , p a +Pb - Pi)t k, 


where Xj ab = z — Vj and Vj = k 2 /(1 — z), cf. Sec. [3j and where ki stands for any final-state momentum, 
including EW particles. The Lorentz transformation, A, is given by 


A (K,K)<i = g»- 


2(K + ky(K + K) v 2 K IJ K V 


(K + K) 2 


K 2 


(A.17) 


The inverted kinematics needed for matching the DIRE shower to NLO computations and for merging with 
higher-order tree-level calculations are given in [3D]. The combined momenta are obtained from Eqs. (A. 1), 


(A.5), (A.7) and (A.13) by swapping momenta and masses before and after emission. In the case of final-state 


emitter with final-state spectator, for example, this amounts to the replacement pu f-A pk and s,y O m 2 -. 


li We changed the definition of the transverse momentum in Eqs. (| A. 1 2|i compared to Ref. [9j, in order to match Eq. (| \. [ 1 [i. 
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B Anomalous dimensions 


The anomalous dimensions of the splitting functions listed in Eqs. (3.4) are computed in this appendix as 

(B.l) 

They are given by 


7 ab(N,K 2 )= [ d ZZ N Pab{z,K 2 ) . 
J 0 


7 n (N, n 2 )=2C F T{N, k 2 )- 
7 gg (N, k 2 ) = 2C f K(N , At 2 )- 


C F (2TV + 3) 

(TV + 1)(TV + 2) ' 

C f (N + 3) 

(TV + 1)(TV + 2) 


lgg{N, K 2 ) = 2C a T(iV, ac 2 ) + 2 C A K(TV, ac 2 ) - 


■7? 


2C a (TV + 3) 2 C A 


(TV + 1)(TV + 2) TV+ 3 


■7 S 


7 qg(N, AC ) — 


where 


TrTV 


+ 


2T fl 


(TV + 1)(TV + 2) TV+ 3 


O p /l 1 3.JV+3 jV+4 . ,„-2f i , 

2T(TV, k 2 ) = 3 2 (:,’ 2 :. , 2 „’ J 2 ^ - In + ' 


2K(TV, At 2 ) = 


(TV + 1)(TV + 2) k 2 

9 P Tl TV+2.AV+4. —2^ 

22^1 (t, -2- ; -ft- ) 

(TV + 2) At 2 


(B.2) 


(B.3) 


By construction, only the soft enhanced terms differ from the DGLAP result. The Altarelli-Parisi splitting 
functions would give T(TV) = —2 H^, with the TVth harmonic number, and K(TV) = 1/TV [241 . 


C Momentum mapping and q T spectra in Drell-Yan type processes 


This section presents a brief phenomenological analysis of the different recoil strategies described in Apps. pOi| 
and A.4 We investigate the impact on the transverse momentum (q F ) spectrum of the Drell-Yan lepton 


pair at the LHC. 

It was pointed out |45l that the non-singlet initial-state parton evolution is generated in dipole-like par- 
ton showers by initial-state splitters with final-state spectator, except for the first branching, which stems 
from an initial-state splitter with initial-state spectator. The kinematics mapping described in App. |A.3| 
then results in the Drell-Yan lepton pair acquiring its entire transverse momentum in the first branching 
Comparing to analytical resummation in impact parameter space m, one is led to conclude that only lead¬ 
ing logarithms can be resummed in this scheme [4fij . However, comparing to resummation in transverse 
momentum space |47j . the modified leading logarithmic structure characteristic for parton showers emerges 


from Eq. (2.10). A kinematics mapping more appropriate for comparison with analytical approaches is given 
by App. E ,4[ Here we restrict ourselves to analyzing its impact on the q F spectrum at the LHC. Similar 


analyses have been performed elsewhere for a variety of other observables using standard dipole-like parton 
showers based on Catani-Seymour subtraction [25J |3Z ■ 

Figure[9 shows DIRE predictions from SHERPA for the two different kinematics schemes described in Apps. pO| 
and |A.4 To highlight the differences in the resummation, we present pure parton shower results. Corre¬ 


spondingly, experimental data are omitted. The global recoil scheme from |A.4| shows a small difference at 
low q F . Similar conclusions were reached in [28> 152], 


4 Note that this argument is incomplete, as it is based on non-singlet parton evolution only. The leading real-emission 
configurations at center-of-mass energies much larger than the di-lepton mass are not the qq —A llg channels, but qg —A llq. 
These processes are enhanced by a large gluon PDF. Similarly, the leading double real-emission configurations are gg —A llqq. 
Both types of configurations contain an initial-initial dipole, which - when radiating additional partons - generates recoil on 
the full final state, including the di-lepton pair. 
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Figure 9: The impact of the kinematics mapping on the qr spectrum of the Drell-Yan lepton pair at LHC 
energies. See Fig.[7]for details on the analysis. We compare the local recoil scheme from App. |A.3| 
with the global recoil scheme from App. A.4 
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